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In this paper, taking in to account Brans-Dick theory, we investigate thermodynamic behavior 
of charged black hole solutions. We study the analogy of the black hole solution with the Van der 
Waals liquid-gas system in the extended phase space by considering the cosmological constant as 
dynamical pressure. We obtain critical values of thermodynamic coordinates and plot P — r+ and 
G — T diagrams to study the phase transition. 


I. INTRODUCTION 

Einstein constructed the theory of general relativity that describes the dynamics of our solar system well enough, 
but it probably does not describe gravity accurately at all scales. The problem that general relativity faced is that, 
it does not accommodate either Mach’s principle or Dirac’s large-number hypothesis. It is unable to describe the 
accelerated expansion of the universe accurately. Herein cosmologists explored various alternative theories of gravity 
[l| . Brans and Dicke were pioneers in studying these alternative theories and they developed another relativistic theory 
known as Brans-Dicke (BD) theory Q. This theory can be regarded as an economic modification of Einstein general 
relativity which describes gravitation in terms of metric as well as a scalar held and it accommodate both Mach’s 
principle and Dirac’s large- number hypothesis. Due to the importance of black holes and gravitational collapse in 
both classical and quantum gravity, authors have investigated various aspects of them in BD theory Q. It has been 
proved that in four dimensions, the stationary and vacuum BD solution is just the Kerr solution with a constant 
scalar held In order to investigate the distinction between the BD theory and Einstein theory Cai and Myung 
proved that the black hole solution in the BD-Maxwell theory in four dimensions is just the Reissner-Nordstrom (RN) 
solution with a trivial scalar held In higher dimensions, however, it would be the RN solution with a non-trivial 
scalar held. This is because the stress energy tensor of Maxwell held is not traceless in higher dimensions and the 
action of Maxwell held is not invariant under the conformal transition. 

On the other hand, thermodynamic properties of the black holes have been fascinating subject for many years. It 
was found out that black holes along all assigned thermodynamic variables also have rich phase structure in complete 
analogy with non-gravitational thermodynamic system similar to van der Waals gas system. With the conception 
of expecting the cosmological constant term to arise from the vacuum expectation value of a quantum held, we can 
assume that it can vary. Hence, we can treat the cosmological constant and its conjugate as dynamical pressure and 
volume of a black hole system respectively . Studying the thermodynamics of black holes in AdS space time has 
exhibit various phase transitions with the same critical behavior as van der Waals model, qualitatively [9| . The paper 
of Hawking was the initial studies on this subject 0. He pointed out there is a thermal radiation (black hole hrst 
order phase transition) for Schwarzschild- AdS black hole space time. Adding charge and/or radiation will result a 
behavior similar to a van der Waals liquid/gas [iMl and the analogy will improve by being in the extended phase 
space where the cosmological constant is interpreted as thermodynamical pressure. 

In this paper, we want to investigate the thermodynamic phase transition of charged black holes in BD theory by 
using the analogy between our system and the van der Waals liquid/gas. 

The outline of our paper is as follows. Section |IT] is devoted to brief review of BD-Maxwell field equations with 
their relation with dilaton gravity by a conformal transformation. In Sec. Ell we obtain charged black hole solutions 
of both dilaton gravity and BD theory. Next, we extend the phase space by considering cosmological constant as 
thermodynamic pressure and calculate critical values, and then we plot diagrams for different cases in Sec. IIVI In 
next section, we give a detailed discussion regarding diagrams, their physical interpretations, and the effects of BD 
parameter. We finish our paper with some closing remarks. 
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II. FIELD EQUATIONS AND CONFORMAL TRANSFORMATIONS 


The action of (n + 1)- dimensional BD-Maxwell theory with a scalar field <f> and a self-interacting potential y($) 
can be written as @ 


Ig = - 


^ - U($) - , 


( 1 ) 


where TZ is the scalar curvature, the factor u is the coupling constant, F^i, = — d^A^ is the electromagnetic 

tensor field and is the electromagnetic potential. Varying the action m with respect to the gravitational fieldo^j/, 
the scalar field $ and the gauge field A^, one can obtain equations of motion with the following explicit forms @ 


r - — 
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2[(n — 1) Lo n] 2 [(n — 1) a; -|- n)] 


(n-l)ci>^^-(n + l)V(d>) 
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( 4 ) 


where and are, respectively, the Einstein tensor and covariant derivative of manifold M. with metric 
Due to the appearance of the second derivatives of scalar field in the right hand side of ([2]), solving the field equations 
©-([I]) directly is a non-trivial task. Using a suitable conformal transformation, one can remove this difficulty. Indeed, 
via the conformal transformation Q the BD-Maxwell theory can be transformed into the Einstein- Maxwell theory 
with a minimally coupled scalar dilaton field. Suitable conformal transformation can be shown as 


^ n - 3 

$ = -ln$, 

4a 


( 5 ) 


where 


= (n — 3)/\/4(n — l)a; -I- 4n. 


( 6 ) 

It is notable that all functions and quantities in Jordan frame {g^v, d’ and F^i,) can be transformed into Einstein 
frame $ and F^^). Applying the mentioned conformal transformation on the BD action ([T]), one can obtain 

action of dilaton gravity 


Ig = - 


IGtt 


jn+l 
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n — 1 


(V<py - U($) - exp - 


4a$ 
(n- I) 




( 7 ) 


where 7Z and V are, respectively, the Ricci scalar and covariant derivative corresponding to the metric g^^y, and U(4>) 


IS 


U($) = (8) 

Regarding (n -|- 1)—dimensional Einstein-Maxwell-dilaton action ([7|), a is an arbitrary constant that governs the 
strength between the dilaton and Maxwell fields. One can obtain the equations of motion by varying this action © 
with respect to g^v, 4* and F^jy 


TZ^i, — 

n — 1 




+ ( F^y^py - 
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V^d> = 


.-IdV 




( 10 ) 


da 




= 0 


( 11 ) 


By assuming the (5^1,, F^,y, $) as solutions of Eqs. (l^- (fTT]l with potential V ($) and comparing Eqs. ©-© with 
Eqs. CT- (fTT]) we find the solutions of Eqs. with potential V ($) can be written as 


Ffj,ut 
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III. CHARGED SOLUTIONS IN (n + 1)- DIMENSIONS 


Our strategy is to construct the solutions of BD theory with n ^ 4 and the quadratic potential 

U($) = 2 A$ 2 . 

Applying the conformal transformation ([5|), the potential U($) becomes a Liouville-type potential 

4a$ 


U($) = 2 Aexp 


n — 1 


(13) 


(14) 


In other words, instead of solving Eqs. ©-O with quadratic potential, we solve Eqs. ©-(mi) with Liouville-type 
potential. Assuming the (n -|- 1)—dimensional metric 


ds^ = -f{r)df -f -f r‘^R^{r)dnl_^, 


(15) 


where dn^_i is the metric of a unit (n — 1)—sphere, and f{r) and R(r) are metric functions. By integrating the 
Maxwell equation (ED, we can obtain the nonzero electric field Ftr as 


Ffr — 


(ri?) 


71 — 1 


exp 
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(16) 


Taking into account the metric (TTSl) with Maxwell field ((T51) . the solutions of © and (nnD are 

f(r) = (»-2) + 2A(a^^l)^ _ m 

’ (a2 + n- 2 )(a 2 _l) (^n - - n) r("-2) 


-b 
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(17) 

(18) 
(19) 


where m is an integration constant which is related to the total mass, c is another arbitrary constant related to the 
scalar field and 7 = a^/(l -b a^). 

Now, we are in a position to obtain the solutions of Eqs. ©-© by using the conformal transformation. Considering 
the following spherically symmetric metric 


ds^ = —U{r)dt^ -b 


dr^ 

V{r) 


+ r^H^{r)dni_, 


( 20 ) 


with Eqs. ©-©, we find that the functions U (r) and V (r) are 
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Using the conformal transformation the electromagnetic field becomes 

„(3-rt)7 


Ftr — 


qc'- 


( 22 ) 


(23) 


j,(n-3)(l-7)-H2 ■ 

As one can see electromagnetic field becomes zero as r —>■ 00 . It is also notable that obtained solutions are just 
the charged solutions of Einstein gravity (Reissner- Nordstrom AdS black hole) as uj —>■ 00 . 
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IV. EXTENDED PHASE SPACE AND P- V CRITICALITY IN BD BLACK HOLES 


Calculations show that the Hawking temperature of a BD black hole on the outer horizon r+ is 


K 

2tt 


47 r V U \ dr 


(24) 


=r+ 


where k is the surface gravity. After some simplifications, we obtain 

^ ^ 2(1 + 0 ^) ^ ( 2 „- 3 )( 7 -l) 

47r(n — 1) \ r+T" 

[7 (n - 3) - n + 2] (a^ + l)^ (n - 2) / c \ 
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(25) 


which is invariant under the conformal transformation because the conformal parameter is regular at the horizon. 
The finite mass and the entropy of the black hole can be obtained by using the Euclidian action Q 


M = 


d"-i)7 / n-I 


IGtt 


1 + 


* - 


(26) 


(27) 


We regard A and its corresponding conjugate quantity as the thermodynamic pressure P = and volume 


V = 


(28) 


respectively 

'dH\ _ fdM\ 

. ) S,Q ) S,Q 

Herein, we are interested in studying the phase transition of this black hole. The equation of state of the black hole 
can be obtain using equation 
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(n — 1 ) / b 


8nr, 
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(29) 


IGtt (q;2 — 1) (q; 2 + n — 2) 

where r_|_ is linear function of the specific volume v in geometric unit Q. 

We can investigate the existence of phase transition and critical behavior of this black hole by plotting and analyzing 
the graphs of P — u and G — T diagrams. One may use the inflection point properties 

'dP'' 
dv 


= 0 , 


rd^p\ 

\ dv"^ Jr 


= 0 , 


to obtain the critical values for the temperature, pressure and volume. Due to the difficulties of solving these equations 
analytically, we use the numerical method to obtain critical values. 

According to first law of black hole thermodynamics and the interpretation of M (total mass of black hole) as PI 
(the black hole enthalpy) [l^ the Gibbs free energy of black hole can be written as 


G = H-TS = M-TS 


(30) 
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FIG. 1: P — r+ (left), G — T (right) diagrams for c=l, n = 4, 5=1 and a; = 1. 

P — r+ diagram, from up to bottom T = 1.8Tc, T = l.STc, T = Tc, T = O.STc and T = 0.5Tc respectively. 
G — T diagram, from up to bottom P = l.SPc, P = Pc and P — 0.5Pc, respectively. 




FIG. 2: P — r+ (left), G — T (right) diagrams for c=l, n = 4, g=l and w = 3. 

P — r+ diagram, from up to bottom T = l.STc, T — l.STc, T = Tc, T = O.STc and T = O.STc respectively. 
G — T diagram, from up to bottom P = l.STc, P — Pc and P — O.STc, respectively. 


The behavior of Gibbs free energy with respect to temperature may be investigated by plotting the graph oi G — T. 
We will see the characteristics swallow-tail behavior which guarantees the existence of the phase transitions. 


V. DISCUSSION ON THE RESULTS OF DIAGRAMS 

Thermodynamical behavior of the system is shown in Figs. [T]0 These critical values represent phase transition 
points, which one can see in G — T and P — r+ diagrams. Studying P — r+ graphs (left panel of Figs. HJIHI) show that 
obtained values are critical points in which phase transition occurs. It can be seen in the graphs that as the coupling 
constant (w) increases, the temperature in which phase transition occurs decreases. Increasing w leads to decreasing 
in the value of Gibbs free energy of phase transition point (see Fig. [3 for more details). The results indicate that with 
larger coupling constant the energy that system needs in order to have phase transition becomes less. It is also evident 
from the effect of the coupling constant on the total finite mass of the black hole that by increasing the coupling 
constant the total finite mass which is interpreted as enthalpy of the system will increase too. It means that in order 
to have phase transition it is expected that the mentioned black hole absorb more mass from surrounding. 

Studying T — r+ diagrams shows that as coupling constant increases, both critical pressure and horizon radius of 
critical point decrease. On the other hand, it is worthwhile to mention that due to the relation between pressure and 
cosmological constant (which is related to the asymptotical curvature of the background), as the coupling constant 
increases, the necessity of having a background with more curvature increases too. 
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FIG. 6: P — r+ (left), G — T (right) diagrams for c=l, n = 6, g=l and cj = 5. 

P — r+ diagram, from up to bottom T = 1.8Tc, T — l.STc, T = Tc, T = O.STc and T = 0.5Tc respectively. 
G — T diagram, from up to bottom P = 1.5Pc, P = Pc and P — 0.5Pc, respectively. 




FIG. 7: P — r+ (left), G — T (right) diagrams for c = 1, n = 4, g = 1. 

P — r+ diagram, for T = Tc, OJ = 1 (solid line), w = 3 (dotted line), cu = 5 (dashed line) and w = 10 (dasheddotted line). 

G — T diagram, for P — 0.5Pc, ca = 1 (solid line), cu = 3 (dotted line), oj = 5 (dashed line) and ca = 10 (dasheddotted line). 




FIG. 8: P — r+ (left), G — T (right) diagrams for c = 1, n = 6, g = 1. 

P — r+ diagram, for T = Tc, w = 1 (solid line), w = 3 (dotted line), ca = 5 (dashed line) and w = 10 (dasheddotted line). 

G — T diagram, for P = 0.5Pc, to = 1 (solid line), ca = 3 (dotted line), cj = 5 (dashed line) and uj = 10 (dasheddotted line). 














FIG. 9: P — r+ (left), G — T (right) diagrams for c = 1, a; = 1, g = 1. 

P — r+ diagram, for T = Tc, n = 4 (solid line), n = 5 (dotted line) and n = 6 (dashed line). 

G — T diagram, for P = 0.5Pc, n = 4 (solid line), n = 5 (dotted line) and n — Q (dashed line). 


Using Fig. IH we can also discuss the effect of dimensionality on critical point and behavior of system. It can be 
seen from the obtained G — T diagrams for different dimensions that as dimensionality increases, the temperature of 
critical point increases too. According to P — r^ diagrams the pressure in which phase transition occurs will increases 
as dimensionality increases. As a consequence of increment in pressure one can see that the cosmological constant 
decreases, due to the relation of P = — ^. Hence the need of having a background with higher value of curvature 
decreases in higher dimensions. It can be inferred that as dimensionality increases the system needs to absorb more 
mass in order to have phase transition. 

At the end one can see the critical horizon radius decreases as coupling constant increases same as critical temper¬ 
ature and pressure. On the other hand for higher dimensions, we have higher values of critical points. The results 
arisen from graphs can also be seen directly through the tables I and II. 


UJ _Tc_7^_ Pc 

I.OOOO 1.5227 0.1869 0.03^ 

3.0000 1.5099 0.1790 0.0371 

5.0000 1.5052 0.1762 0.0366 

10.0000 1.5010 0.1736 0.0361 

Table I: critical quantities for <7 = 1 and n = 4. 


UJ _Tc_7^_ Pc 

1.0000 1.2859 0.7460 0.3o1a 

3.0000 1.2456 0.5839 0.2627 

5.0000 1.2325 0.5517 0.2434 

10.0000 1.2210 0.5130 0.2318 

Table II: critical quantities for <7=1 and n = 6. 


VI. CONCLUSIONS 

In this paper, we considered the BD theory in the presence of electromagnetic field and studied its phase structure. 
We extended the phase space by considering cosmological constant as thermodynamical pressure and its conjugate 
variable as volume and regarded the interpretation of total mass of black hole as the enthalpy of the system. 

Studying calculated critical values through two different types of phase diagrams resulted into phase transition 
taking place in the critical values. Studying P — r_|_ and G — T diagrams exhibits similar behavior near critical points 
to their corresponding diagrams in Van der Waals liquid/gas. 
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The results indicated that for large values of coupling constant the system needs less energy (mass) absorption to 
have phase transition, due to the fact that as coupling constant increases the critical temperature decreases. On the 
other hand, studying the effects of dimensionality showed that for higher dimensional black holes, phase transition 
take places in higher temperature and lower Gibbs free energy. 

Finally considering BD theory with various models of nonlinear electrodynamics, it would be interesting to analyze 
the effects of nonlinearity on extended phase space thermodynamics and P —V criticality of black hole solutions. We 
left these issues for the forthcoming work. 
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